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Question 1 [10 Marks] 


(a) 


A vector is input in MATLAB via the command: 
aS 246d 


Write down the expected output from the following subsequent commands: 


Consider the following function, func.m: 


function [f] = func(x,y) 
f=x"2+y~2; 


end 


What would be the associated output from the MATLAB command z=func (1,2)? 


Assuming that x is a vector of real numbers, what does the following MATLAB 


function do? 


function m = func2(x) 
n=length(x) ; 
m=x (1) ; 
for 11=2:n 

if x(ii)>m 

m=x (ii); 

end 

end 


end 
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Question 2 [10 Marks] 


(a) Write a MATLAB function to determine the average of all the elements of a matrix 


mat without using MATLAB’s built-in mean command. 


(b) Write a MATLAB script to evaluate y = tan at 1000 equally spaced values of x 
over the interval —~ ao a) and then graph the function as a red dashed line. 
Add code to your script to label both the x- and y-axes of the graph, and to add a 
title. 
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Question 3 [10 Marks] 


(a) Convert the following numbers from binary to decimal form: 

(i) 101, 

(ii) 0.101, 
(iii) 101.101, 
(iv) 110.0011, 

(b) Convert the following numbers from decimal to binary: 

(i) 73 

(ii) 0.75 

(iii) 73.75 
(c) Re-write the binary representation of 73.75 as a double precision floating point num- 

ber (in binary) according to IEEE standard 754. 
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Question 4 [5 Marks] 


(a) 


The general Taylor series expansions for a function, f(x), expanded about a point 


X = Xp, and evaluated at the points x,,, and x,_, can be written as: 


h? h3 ht : 
J anaa) = Fy) a7 by (ie) ae af (en) a gif (en) he af an) 7 O(h”) 


and 
he he eer 
Fina) = flim) =AS On) + hon) = SF" en) + FLOMan) +O) 
respectively, where the points x = ..., %p_-1, Un, Ln41,... are equally spaced with sep- 


aration h = %pj41 — Lp, and O(h?) indicates terms that are of the order of h°® and 
above. Use the above Taylor series to obtain approximations for f’(x,) and f” (an) 


that have associated errors of O(h?). 


Another standard form for the Taylor series of a suitable function f(x) about a point 


"a "a (n) a 
Flo) = f(a) + FO — 0) + AO (a — ay? +. + Fo — a)" + yal), 
where f(y.) 
— Ln, —a n+1 
Rn+i(Z) (n+ 1)! (x ) ’ 


and fl, is Some number between x and a. Use the above form of the Taylor series to 


determine the first three non-zero terms of the Taylor series expansion of sin x about 
sin x 


x = 0. Use the terms of your series to then find an approximation of when «x is 


close to zero. 
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Question 5 [5 Marks] 
The following is an incomplete section of MATLAB script to apply the bisection algo- 


rithm to find a zero of the function contained in a file func3.m within the interval [2, 7,]. 
Complete the indicated lines of code so that the zero is determined correct to 6 decimal 


places. 


clear 
xl=input(’Enter left end point of interval containing solution: ’); 


xr=input(’Enter right end point of interval containing solution: ’); 


nmax=30 ; 
f1=func3 (x1) ; 
for ii=1:nmax 
xnew=... {Incomplete line 1 


fnew=... %Incomplete line 2 


if fl*fnew>0 

xl=xnew; 

fl=... %Incomplete line 3 
else 

xr=... %Incomplete line 4 


end 


if abs(xr-xl)<... %Incomplete line 5 
break 
end 


end 
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Question 6 [10 Marks] 


Consider the function g(x) = cos? x. 


(a) 


Sketch a simple graph of y = g(x) over the interval 0 < x < 27, along with a graph of 
y = x. Based on your graph, identify an approximate solution to x = g(x) that could 
be used as the initial input, xo, into a fixed point iteration algorithm to determine a 


more accurate solution. 


What is the general condition, in terms of g/(x), for fixed point iteration to converge 
to a solution of x = g(x). Given that the root of x = cos? x is approximately 0.642, do 
you expect a fixed point iteration algorithm to converge to the solution of x = cos? 2, 


or not? 


Re-write the equation x = g(x) in the form f(x) = 0 (explicitly including the cos? x 
term). Then, write down the Newton-Raphson iteration formula for solving this 


specific problem. 


The following is an incomplete MATLAB script to determine the solution to the 
specific equation x = cos? x using Newton-Raphson iteration. Complete the indicated 
lines of code so that the solution is determined, correct to 7 decimal places, and with 


a maximum of 20 iterations. 


xold=0.75; 
for ii=... %Incomplete line 1 
xnew=... 4Incomplete line 2 


if ...<5e-8 %Incomplete line 3 
break 

end 

...=xnew; ZIncomplete line 4 


end 
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Question 7 [10 Marks] 
Consider the following data: 


x |-1 | -0.96 | -0.86 | -0.79 | 0.22 | 0.50 | 0.93 
y | -1 | -0.151 | 0.894 | 0.986 | 0.895 | 0.500 | -0.306 


(a) Draw a graph that illustrates the approximate positions of the given data points in 


the xy-plane. 
(b) What is the degree of the polynomial that passes exactly through all the data points? 


(c) If you were to estimate y(0.75) using a linear approximation, which pair of points 


would be the most appropriate to use? 


(d) The MATLAB command, p = polyfit(x,y,n) returns a vector p of coefficients of 
the least squares polynomial of degree n interpolating the points (x(k) ,y(k)) for 
k=1, 2, ..., n+1. The related command polyval(p,z) returns the value of the 
interpolating polynomial with coefficients p at the point z. With the aid of these two 
commands, write a short MATLAB script to estimate y(0.75) using a quadratic fit 


to the data given above. 


(e) Explain how a cubic spline fit to a set of data is constructed. Based on your graph 
in part (a) of this question, do you think that a cubic spline approximation would be 


an appropriate method of curve fitting for the data given in the table above? 


Please remember: This examination paper MUST BE HANDED IN. Failure to do 


so may result in the cancellation of all marks for this examination. Writing your 
name and number on the front will help us confirm that your paper has been 


returned. 


